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We study the electronic properties of strongly spin-orbit coupled electrons on the elastic pyrochlore 
lattice. Akin to the Peierls transition in one-dimensional systems, the coupling of the lattice to 
the electronic degrees of freedom can stabilize a spontaneous deformation of the crystal. This 
deformation corresponds to a breathing mode, which breaks the inversion symmetry. We find 
that for intermediate values of the staggered strain, the inversion-symmetry broken phase realizes a 
topological Weyl semimetal. In the temperature-elasticity phase diagram, the Weyl semimetal shows 
a reentrant phase behavior: it can be reached from a symmetric phase realized both at higher and 
at lower temperatures. The symmetric phase is a Dirac semimetal, which is protected by the non- 
symmorphic space group of the pyrochlore lattice. Beyond a critical value of the staggered strain, the 
symmetry-broken phase is a fully gapped trivial insulator. The surface states of the Weyl semimetal 
form open Fermi arcs and we observe that their connectivity depends on the termination of the 
crystal. In particular, for the {111} films, the semiclassical closed electronic orbits of the surface 
states in a magnetic field cross the bulk either twice, four, six or twelve times. We demonstrate 
how one can tune the number of bulk crossings through a Lifshitz-like transition of the Fermi arcs, 
which we call Weyl-Lifshitz transition, by applying a surface potential. Our results offer a route to 
a topological Weyl semimetal in nonmagnetic materials and might be relevant for pyrochlore oxides 
with heavy transition-metal ions such as alloys of iridates. 


I. INTRODUCTION 

Pyrochlore oxides experience growing interest as po¬ 
tential hosts of novel electronic phases. Due to the pres¬ 
ence of both strong spin-orbit coupling and electronic cor¬ 
relation effects, pyrochlore oxides with heavy transition- 
metal ions have been identified as a playground for topo¬ 
logical phases [1]. Some of the notable theoretical predic¬ 
tions include exotic spin liquids [2, 3], topological insu¬ 
lators [4-6], topological crystalline insulators [7], topo¬ 
logical semimetals [8, 9] and unusual forms of mag¬ 
netism [10, 11]. Besides these intrinsic three dimensional 
phases, heterostructuring of pyrochlore oxides offers the 
possibility to access two-dimensional topological phases, 
such as the quantum spin Hall or Chern insulator [12-14]. 

In this work, we add an additional ingredient to the 
physics of pyrochlore oxides: the coupling of the elec¬ 
tronic degrees of freedom to the lattice. Our starting 
point is a general Hamiltonian on the pyrochlore lat¬ 
tice, which, as a consequence of the non-symmorphic 
space group, realizes a Dirac semimetal with a discrete 
set of four-fold degenerated Fermi points. We find that 
for an intermediate stiffness of the crystal, an inversion- 
symmetry breaking staggered strain can spontaneously 
develop. The symmetry-broken phase realizes a topolog¬ 
ical Weyl semimetal, whose Fermi surface consists of a 
discrete set of doubly degenerated Fermi points with lin¬ 
early touching valence and conduction bands. If the stag¬ 
gered strain reaches a critical value, the system turns into 
a fully gapped trivial insulator. Our model shows a rich 
phase structure including a reentrant Weyl semimetal, 
which can be reached from the symmetric phase present 
at both low and high temperatures. 

Weyl semimetals have originally been proposed by 
Wan et. al. [15] in magnetic pyrochlore iridates, such as 


Y 2 lr 2 07 . They have been the subject of several reviews 
[16-18]. Recent investigations of Weyl semimetals fo¬ 
cused on their realizations in interacting models [8, 9, 15] 
and in normal insulator-topological insulator multilay¬ 
ers [19-21]. Recently, Weng et. al. [22] proposed that 
certain non-centrosymmetric transition metal monophos¬ 
phides might be nonmagnetic Weyl semimetals. 

The exciting properties of Weyl semimetals are diverse: 
First, they host unusual surface states which form dis¬ 
joint open Fermi arcs rather than closed Fermi lines. The 
Fermi arcs connect the projections of the Weyl nodes onto 
the surface Brillouin zone and Fermi arcs on opposite sur¬ 
faces of the sample are linked through bulk states. These 
properties lead to characteristic signatures in quantum 
oscillations experiments in a magnetic field [23]. Second, 
also the bulk properties of Weyl semimetals are unusual 
due to the chiral anomaly, which states that the elec¬ 
tric charge carried by electrons of a given chirality is not 
conserved in the simultaneous presence of electric and 
magnetic fields. In contrast, the Nielsen-Ninomiya dou¬ 
bling theorem guarantees that the net current carried 
by electrons of both chiralities is conserved [24]. The 
chiral anomaly implies unusual transport properties, re¬ 
viewed by Hosur and Qi [17] and by Ramamurthy and 
Hughes [18], which include the axion response [25], the 
semiquantized anomalous Hall effect [26, 27], and the 
(still controversial) chiral magnetic effect [25, 28, 29]. 
A direct way to probe the chiral anomaly in topologi¬ 
cal semimetals has also been proposed using a non-local 
transport experiment [30]. 

The Weyl and Dirac semimetals arise naturally as the 
intermediate phase in a topological insulator - normal in¬ 
sulator phase transition [31-35]. More generally, Weyl 
semimetal can be obtained by starting with a three- 
dimensional Dirac semimetal and then breaking inver- 
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sion or time-reversal symmetry [19-21, 36]. Recent in¬ 
vestigation led to the discovery that NaaBi and CdsAs2 
are Dirac semimetals [37-41]. Young et. al. showed 
that Dirac semimetal can arise in non-interacting mod¬ 
els as a symmetry protected phase [42] due to the non- 
symmorphic elements of the crystal space group. Non- 
symmorphicity has recently been found to have further 
effects on the electronic structure, for example, it pro¬ 
hibits band insulators at certain integer fillings [43] or can 
protect a novel Z2 topological crystalline insulators [44]. 

Our paper is organized as follows. In Sec. II, we intro¬ 
duce a general tight-binding model of spin-orbit coupled 
electrons on an elastic pyrochlore lattice. In Sec. Ill we 
discuss how the deformation of the crystal affects the 
band structure, and we show phase diagrams in the ab¬ 
sence and presence of a symmetry-breaking staggered 
stress field. Due to the importance of the space group 
in realizing these intriguing electronic phases, we present 
in Sec. IV a detailed group theoretical discussion of the 
spectrum, which particularly highlights the role of non- 
symmorphicity. In Sec. V we discuss the surface states of 
{111} and {111} oriented films. We show that the surface 
states depend on the termination of the lattice and intro¬ 
duce the concept of the Weyl-Lifshitz transition, which 
is characterized by a change in the connectivity of the 
Fermi arcs. 


II. MODEL 

A. Spin orbit-coupled electrons on the pyrochlore 
lattice 

We study a system of non-interacting electrons on the 
lattice of corner-sharing tetrahedra called the pyrochlore 
lattice. These tetrahedra form a bipartite lattice and, 
thus, can be labelled as even and odd [illustrated in 
bright blue and dark orange in Fig. 1(d)] in such a way 
that every even tetrahedron touches only odd ones and 
vice versa. The pyrochlore lattice is a face-centred cu¬ 
bic (FCC) lattice belonging to a non-symmorphic space 
group #227 (Fd3m). Its Brillouin zone has the form of 
a truncated octahedron illustrated in Fig. 1(b). 

We consider the following Hamiltonian to describe the 
dynamics of the electrons on the pyrochlore lattice [9] , 

H 0 = ^2 cj (t 1 + it 2 dij • r) Cj + 

(hj) 

+ 'y ^ cj [l| + i (t 2 'R'ij + ’ T ] c j • (1) 

«hj» 

The electron annihilation operator at site i is given by 
c i — (c it ,c 4 ) T , where a refers to a (pseudo-)spin- 

1/2 degree of freedom. The meaning of the symbols is 
as follows: The first sum runs over the nearest neighbor 
(NN) and the second sum over the next-nearest neighbor 
(NNN) bonds. The Pauli matrices r = (t x ,t v ,t z ) de¬ 
scribe a global basis in spin space, and the terms contain¬ 


ing them represent spin-orbit coupling. The real space 
vectors appearing in the NN part of Ho are defined via 


dij = 2aij x bij (2a) 

— 2 T hj) x c (2b) 

bij — bj bi (2d) 

= | (hi + b 2 + 63 + 64) (2e) 

and those in the NNN part are defined as 

'R'ij — bik x h/cji (3b) 

'Dij — difc x dfcj (3c) 


where k is the common NN of the NNN sites i and j. We 
further use position vectors bi pointing to the site i of 
the tetrahedron 


61 

= a(0,0,0) 

(4a) 

*>2 

= a(0,1,1) 

(4b) 


= a(l, 0,1) 

(4c) 

&4 

= a( 1,1,0) 

(4d) 

and position vector x c pointing to the centre 

i of the tetra- 


hedron 

1 4 

x c = b i (5) 

^ i=l 

Examples of bij, aij and d^ are illustrated in Fig. 1(a). 
We denote the length of the cube circumscribed to the 
tetrahedra [Fig. 1(a)] as a so that \bij\ = a\[2. The edge 
length of the FCC unit cell is 4a. 

Hamiltonian (1) is the most general single-orbital 
Hamiltonian with spin-orbit coupling and up to NNN 
terms respecting the full symmetry of the pyrochlore lat¬ 
tice [6]. Such a situation arises for the five 5d-electrons 
coming from each Ir 4+ site of the pyrochlore oxides. The 
crystal field of the neighbouring atoms splits the 5 d or¬ 
bitals into the six t 2g states and the four e g states, the 
latter being higher in energy. The on-site spin-orbit cou¬ 
pling further splits the six degenerate t 2g states into an 
effective pseudospin j e ff = 1/2 doublet and an effective 
Jeff = 3/2 quadruplet, the latter being lower in energy. 
The five 5d-electrons of Ir 4+ sites fill the j e ff = 3/2 
quadruplet and half-fill the j e q = 1/2 doublet [3]. Hamil¬ 
tonian (1) can then be viewed as an effective model for 
electrons residing in the j e ff =1/2 orbitals. In this way 
one can relate the parameters ti :2 and t[ 2 3 to particular 
orbital overlaps as is thoroughly explained in Ref. [9]. 

A generic spectrum of Hamiltonian (1) is plotted in 
Fig. 2(a). Note that the site fillings n = 1/4 and n = 3/4 
correspond to a semimetallic phase. In pyrochlore iri- 
dates, such a commensurate filling may be realized by 
considering alloys of the form A2-^B cc Ir207 where A and 
B are nonmagnetic but have different oxidation states 
[45-47]. For example, we expect that Yi.sCao.51^07 
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FIG. 1. (Color online) (a) Illustration of vectors x c , a,ij , 
bij , and dij used in Hamiltonians (1) and (9). (b) Shape of 
the Brillouin zone (BZ) of a face-centred cubic (FCC) lat¬ 
tice. This corresponds to both the pyrochlore lattice and 
the breathing pyrochlore lattice, (c) Definition of the high- 
symmetry points of the FCC lattice within one eighth of the 
Brillouin zone and the path rXWLTKX in the momentum 
space used for plotting energy spectra in Fig. 2. The U point 
is equivalent to the K point, (d) Pyrochlore lattice with differ¬ 
ently colored even (blue/bright) and odd (orange/dark) tetra- 
hedra. (e) The Z-broken (breathing) pyrochlore lattice con¬ 
sidered in the elastic model. 


[45] realizes the site filling n = 1/4 and the (hypotheti¬ 
cal) compound Bii.sSeo.sI^Oy site filling n = 3/4 in our 
model. To be concrete, we will consider 

n = 3/4 (6) 


throughout the paper. 


B. Elastic lattice 

Electron-phonon coupling can lead to softening of cer¬ 
tain phonon modes and to distortion of the lattice akin 
to the Peierls transition. This is also similar to certain 
valence bond solids arising due to interactions [48]. The 
leading instabilities can be found by investigating the 
Lindhard function. We argue that the leading lattice in¬ 
stability of model (1) occurs at momentum q = 0 . To 
illustrate this, we ignore the two occupied bands far from 
the chemical potential and consider only the conduction 
and the valence bands. Defining £ a ,fe = £ a ,k — /i where 
a stands for the band index, Fig. 2(a) indicates approxi¬ 
mate electron-hole symmetry 


(Ccon ,k AO ~ (£val,fe AO (7) 


temperature satisfies 


0<-x(q) = ~7, X X 


O 

fee BZ ol, oc' 

1 2 

12 ik-o 4 

feeBZ q 


f{£,a,k-q ) ~ /(&*',fc) 

a,k — q Ca',fe 


£fe 


- O 


12 


feeBZ 


1 


2_W_ 


2 £k-q ^fk) 


= -x(0) (8) 


where we used the arithmetic-harmonic mean inequality. 
The equality sign applies only if q = 0 or if both bands 
are perfectly fiat. The result means that the Lindhard 
function has a peak at q = 0 which corresponds to the 
leading instability. We expect this peak to be preserved 
for Hamiltonian (1) where the electron-hole symmetry is 
approximately valid. 

Inspired by this observation, we consider the most sym¬ 
metric q = 0 phonon mode which corresponds to the si¬ 
multaneous expansion of the even and shrinking of the 
odd tetrahedra as illustrated in Fig. 1(e). We call it the 
breathing mode of the pyrochlore lattice. Similar breath¬ 
ing pyrochlore lattice has been observed in certain A-site 
ordered spinel oxides [49, 50]. 

We treat the breathing mode classically and refer to its 
amplitude as staggered strain 5. We model its effect on 
the electron Hamiltonian by multiplying the NN terms by 
a factor (1 + (— )S) for the short (long) bonds. We ignore 
the higher order influence on both the NN and the NNN 
terms as these are assumed to have a quantitative but 
not qualitative effect on the phase diagram of the model. 
The elastic energy of the lattice deformation is set to 
be proportional to the square of the amplitude 5. The 
complete Hamiltonian then reads 

Us = ^YQS 2 + ^ (1 ± <J) c\ (G + it 2 d i j • r) c j 
(hj) 

+ 4 [f| + i (t 2 'R'ij + t^'Dij) • t] Cj (9) 

«bi» 

where Y is elasticity of the lattice analogous to the 
Young modulus, and £2 is the volume of the sample. 
The equilibrium value of S is determined by minimizing 
the energy of the electron-lattice system which, by the 
Hellmann-Feynman theorem, corresponds to solving the 
self-consistency equation 


(5 


X T 4 (ii + it 2 dij ■ t) Cj |^) 

<iJ> 


( 10 ) 


where |4G) is the ground state of the electron Hamilto¬ 
nian for a staggered strain amplitude S. 

Hamiltonians (1) and (9) have many free parameters. 
To simplify the situation, we reduce the parameter space 
by setting 


In the case of a perfect electron-hole symmetry £ C on,fe = 
—£vai,fe> the static (cj = 0) Lindhard function at zero 
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for the relative strength of the spin orbit coupling and 


P ' h t 2 


(lib) 


for the relative strength of the NNN terms. We further 
define a variable 


s = sign(ti) = ±1. (He) 

These substitutions modify Hamiltonian (9) to 

Hs = \ Y ™ 2 + s { X ( x ± v 4 ( x + iRd a • T ) c i 

(h3) 

+ P X 4 t 1 + iR + V ij) • T ] C j } ( 12 ) 

«i.J» 

where all energies are now in units of \ti\. We will work 
with 


R = -0.4 (13a) 

p = —0.1 (13b) 

5 = +1 (13c) 

unless otherwise stated. This choice approximately cor¬ 
responds to the parameters used in Ref. [9]. 

It may happen that the energy of the Hamiltonian (12) 
is minimized for a non-zero value of 5. This means that 
the cost of the deformation is compensated for by the 
electron energies, thus making the lattice deformation 
energetically favorable. Such a transition decreases the 
symmetry of the lattice. In particular, it breaks the in¬ 
version symmetry X. The Z-broken lattice still has the 
FCC Bravais lattice and an unchanged Brillouin zone but 
it belongs to a symmorphic #216 (F43m) space group. 


III. PHASE DIAGRAMS 
A. Evolution of the band structure 

We start by describing the evolution of the spectrum 
of Hamiltonian (12) as we tune the staggered strain am¬ 
plitude S. The following observations are based on a di¬ 
rect numerical diagonalization of the Hamiltonian and 
are summarized in Fig. 2. In Sec. IV, we will present 
detailed group theoretical arguments. 

In the absence of a staggered strain, S = 0 [Fig. 2(a)], 
we find that all bands are doubly degenerate as a conse¬ 
quence of the simultaneous presence of time-reversal T 
and inversion symmetry X. At each X point [see Fig. 1(c) 
for the definition of the high-symmetry points], we find 
two energetically separated Dirac nodes where four bands 
reach the same energy and disperse linearly in all direc¬ 
tions. This is similar to the spectrum of graphene but 
in three rather than just two dimensions. Since the in¬ 
equivalent X points are related by crystal symmetries, 
all higher (and all lower) Dirac nodes are realized at the 



FIG. 2. (Color online) The spectrum of Hamiltonian (12) 
for parameters (13) and (a) 5 = 0, which corresponds to the 
Dirac semimetal (DSM), (b) 5 = 0.08, which corresponds to 
the Weyl semimetal (WSM) with 12 distinct Weyl nodes, and 
(c) 5 = 0.13, which corresponds to a trivial band insulator 
(INS). Left figures show the spectra in the three cases along 
the path in the BZ indicated in Fig. 1(c), right figures show 
the spectra over the square face of the BZ. The arrows indicate 
position of one of the Weyl nodes for the WSM phase, and 
the place where the Weyl nodes have annihilated for the INS 
phase. 


same energy. It is therefore possible to tune the chemical 
potential to this value, which corresponds to a site filling 
of n = 3/4 (six electrons per unit cell) for the upper and 
n = 1/4 (two electrons per unit cell) for the lower Dirac 
nodes. Such a system has a Fermi surface consisting of 
a set of k -points and is usually referred to as a Dirac 
semimetal (DSM) [31, 37-42]. 

A non-vanishing staggered strain, 5 # 0 [Fig. 2(b)], 
leads to a splitting of each Dirac node into four Weyl 
nodes. These are points where only two rather than four 
bands touch each other, and they disperse linearly in all 
directions around this point. As we increase the value 
of 5, the Weyl nodes move along the XW lines. The 
upper 12 Weyl nodes of the model are mutually related by 
crystal symmetries and time-reversal. As a consequence, 
they are all realized at the same energy and the chemical 
potential resumes to be tuned to them if either n = 1/4 
or n = 3/4. Such a phase is called the Weyl semimetal 
(WSM) [15-18]. 

In the vicinity of each Weyl node, we can approximate 
the electron Hamiltonian as 


H{k) = (E 0 + Uiki)t + VijkiCTj 


(14) 
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FIG. 3. (Color online) Phase diagrams of Hamiltonian (12) 
for p = —0.1 and s = +1 in the R — Y plane for (a) T — 0, 
and (b) T = 0.5. The red (vertical) arrow in diagram (b) 
corresponds to the same set of parameters as the (horizontal) 
red arrow in Fig. 4(a). At T = 0 (T = 0.5), the transition 
from the DSM phase to the WSM phase is first (second) order 
for the shown parameters. 

where a = (cr x ,cry,<j z ) are the Pauli matrices, det [vij] ^ 
0 and vector u describes the tilt of the dispersion cone. 
The quantity sign(det[u]) = ±1 associated with each 
Weyl node is called chirality. It gives rise to the topolog¬ 
ical nature of the WSM, which manifests itself, e.g., by 
the appearance of the Fermi arcs in the surface Brillouine 
zone [15] and by the chiral anomaly [17]. The robust¬ 
ness of the WSM against all local perturbations follows 
from the fact that all Pauli matrices are used in the ef¬ 
fective low energy Hamiltonian (14) [15]. Gapping the 
spectrum requires enlarging of the effective low-energy 
Hilbert space and can be achieved either by scattering 
between different Weyl nodes (which breaks translational 
symmetry) or by forming a superconducting state (which 
requires breaking the global U( 1) symmetry) [17]. 

Time-reversal symmetry implies that if there is a Weyl 
node at momentum fc, then there is a Weyl node of the 
same chirality at —k. On the other hand, Weyl points 
at momenta related by mirror symmetry carry opposite 
chiralities. These requirements allow us to characterize 


all twelve Weyl nodes. We find 6 Weyl nodes with pos¬ 
itive and six with negative chirality, in accordance with 
the Nielsen-Ninomiya doubling theorem [24]. 

At a critical value of the staggered strain S = S c , pairs 
of Weyl nodes with opposite chirality meet at the W 
points where they annihilate. The spectrum disperses 
quadratically along the XW lines, which correspond to 
the direction of motion of the Weyl nodes [32]. Such 
a quadratic band touching point does not carry chiral¬ 
ity [35]. Finally, for S > S c , we find a gapped phase, 
which is a topologically trivial band insulator (INS). For 
the chosen parameters (13) and filling factor (6), the crit¬ 
ical staggered strain is 

5 C « 0.1112. (15) 


B. Phase diagrams for elastic lattices 

To find the phase diagrams at filling n = 3/4, we 
numerically solve the self-consistency equation ( 10 ) for 
varying hopping parameters and elasticities of the lat¬ 
tice. In general, we expect that the deformation of a stiff 
lattice (i.e. with a large value of Y) is energetically too 
costly, so that the system will remain in the X-preserving 
DSM phase. Decreasing the elasticity should allow for a 
transition to an X-broken state that can be either a WSM 
or an INS. 

We first study the role of the relative spin-orbit cou¬ 
pling i7, see Figs. 3(a) and (b). Throughout this sub¬ 
section, we fix p = —0.1 and s = +1 and we vary R 
in the range (—0.774, +0.364) for which the undistorted 
system realizes a DSM phase. For values of R outside 
of this interval, the valence bands at T rise above the 
Dirac nodes at X, leading to a metallic state with a hole¬ 
like Fermi pocket at the T point and three electron-like 
Fermi pockets at the inequivalent X points. The value 
R = 0 corresponds to no spin-orbit coupling and hence 
spin-independent hopping amplitudes. In this case, the 
spectrum does not exhibit Dirac nodes but line nodes 
along the XW lines. This line degeneracy is gapped out 
by any S ^ 0. The transition from R < 0 to R > 0 also 
changes the degeneracies at the T point (bottom-to-top) 
from 2-4-2 to 4-2 -2. 

At zero temperatures, see Fig. 3(a), the transition from 
S = 0 to S 7 ^ 0 is always first order. Furthermore, if 
R G (—0.35, 0.28)\{0}, a narrow WSM phase is found by 
varying the elasticity. (Our numerical accuracy does not 
allow a definite conclusion on the presence of the WSM 
phase for \R\ < 0.02). For values of R outside of the 
mentioned interval, we find a direct first-order transition 
between the DSM and the INS phase, i.e., the staggered 
strain directly jumps to a value 6 > S c . 

In an equivalent calculation at T = 0.5, the transition 
is found to be second order for all values of R. By re¬ 
ducing Y, S smoothly increases from 0 in the DSM to a 
value 6 > 5 C in the INS, leading to a larger WSM region 
than at zero temperature. 
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FIG. 4. (Color online) (a) Phase diagram of Hamiltonian (12) for R = —0.4, p = —0.1, s = +1 in the R — T plane. The 
WSM phase is most robust at intermediate temperatures. The thick (green) line indicates the first-order transition between the 
symmetric DSM and the symmetry-broken WSM and INS phases and ends at a tricritical point TP. At higher temperatures, 
the transition is always second order. The horizontal (red) arrow corresponds to the same set of parameters as the vertical 
(red) arrow in Fig. 3(b). The blue rectangle corresponds to the range of temperatures and elasticities that are used in the three- 
dimensional diagram, (b) Phase diagram for the same parameters with included coupling to a symmetry-breaking staggered 
strain YS o. In the presence of such a field, the distinction between the DSM and the WSM phases ceases to have a meaning. 
Only the structural first-order transition (the dark green sheet in the back) and the transition between the WSM and INS 
phases (the pale blue sheet in the front) are therefore present. Where the two transitions occur at the same parameters, the dark 
red sheet is plotted. (c,d) Horizontal cuts of the three-dimensional phase diagram in (b) corresponding to fixed temperature, 
and (e,f) vertical cuts corresponding to a fixed staggered strain. In these four diagrams, the dotted green lines correspond to 
the first-order structural transition, the dashed blue lines to the WSM-to-INS transition, and the solid red lines is drawn where 
the two transitions coincide. 


We further study the phase diagram in the Y-T plane 
for fixed R = —0.4 [Fig. 4(a)]. The diagram demonstrates 
that the WSM phase is indeed most robust at intermedi¬ 
ate temperatures. We further observe a reentrant phase 
behavior: Starting in the symmetric Dirac semimetal at 
high temperatures, the system spontaneously breaks the 
inversion symmetry upon cooling. But upon further cool¬ 
ing, it returns to the symmetric phase. For the cho¬ 
sen parameters, we observe this behavior in the range 
Ya 3 G (0.513,0.594). The reentrance can be traced back 
to the peculiar form of the density-of-states (DOS) in 
the symmetric DSM, which has a local minimum with 
a vanishing DOS at the Fermi energy. Thermal broad¬ 
ening then enhances the effective DOS at the chemical 
potential, rendering the system more susceptible to a 
symmetry-breaking transition at elevated temperatures 
than at low temperatures. A similar reentrant phase 
behavior has also been observed in models for metallic 
metamagnetic systems [51]. 

At low temperatures, the transition between the sym¬ 
metric and the symmetry-broken phase is first order. The 
first-order transition line obeys the Clausius-Clapeyron 
relation, and has to approach zero temperature perpen¬ 
dicularly. The first-order transition might also lead to a 
hysteretic behavior with temperature. Slightly below the 
tip of the belly-shaped WSM phase, there is a tricritical 
point (TP). For temperatures above the tricritical point, 


the transition between the symmetric and the symmetry- 
broken phase is second order. 


C. Phase diagram in a symmetry-breaking field 

To further explore the interesting structure of the 
phase diagram, we study the effect of a symmetry¬ 
breaking field, which is conjugate to the staggered strain 
(i.e. a staggered stress). We incorporate the symmetry¬ 
breaking field in the elastic lattice model (12) by modi¬ 
fying the elastic energy term as follows, 

Yns 2 ^ Wn( 5 -s 0 )\ ( 16 ) 

where YS o parametrizes the staggered stress. The result¬ 
ing three-dimensional phase diagram in the (Y,So,T)~ 
space is shown in Fig. 4(b). Due to the complexity of 
the three-dimensional diagram in Fig. 4(b), we show in 
Figs. 4(c)-(f) several slices at either a fixed temperature 
or a fixed staggered stress. 

A non-vanishing explicitly breaks the inversion sym¬ 
metry. Therefore, the DSM changes into a WSM phase 
with a small separation of the Weyl nodes proportional 
to S 0 . This means that the second-order transition line 
separating DSM from WSM at = 0 ceases to ex¬ 
ist in the presence of a staggered stress. However, the 
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first-order line below the triciritical point survives also 
for 5 0 > 0, forming a sheet of first-order transitions (a 
so-called Griffiths wing), which extends up to a finite 
value of Jo- The Griffiths wing either signals a structural 
transition within the WSM phase [dark green sheet in 
Fig. 4(b) in the online version] or between the WSM and 
the INS phases [dark red sheet in Fig. 4(b) in the online 
version]. Note that the Griffiths wing takes an unusual 
form with a “belly” at finite temperatures: the end point 
at T = 0.1 occurs at Jo ~ 0.0030 while at T = 0 it oc¬ 
curs at ~ 0.0016, see Figs. 4(c) and (d). Finally, we 
note that the boundary between the WSM phase and the 
INS phase remains well-defined for all values of Jo and T 
[light blue sheet in Fig. 4(b) in the online version]. 


IV. GROUP THEORETICAL ANALYSIS 

A. The general strategy 

In this section, we will demonstrate how the symme¬ 
try of the pyrochlore lattice inevitably leads to the Dirac 
node at the X point of the Brillouin zone, and why this 
Dirac node has to split into four Weyl nodes upon break¬ 
ing the inversion symmetry. The pyrochlore lattice be¬ 
longs to the same space group as the diamond lattice 
(#227, Fd3m) so the same reasoning as that of Ref. [42] 
applies. Further structures belonging to this space group 
are /3-cristobalite and spinel oxides. 

Our main tool in determining the spectrum degenera¬ 
cies at a given k- point is the following relation between 
representations of point-symmetry operations in the vec¬ 
tor space spanned by the Bloch wave functions at k: 

D k {Ri)D k {Rj) = exp (~ig i • tj)D k (Ri o R 3 ). (17) 

Here, Ri is a point-symmetry operation that maps the 
considered fc-point onto itself modulo a reciprocal lattice 
vector g i = (R^k) — k , and t{ is the non-symmorphic 
shift associated with the point operation Ri. The set 

_ k 

of all such /c-preserving operations Ri forms a group G 
called the little co-group of k. The function 

6{Ri, Rj ) = exp (-igi • tj) (18) 

is called the factor system of the representation and it is 
completely fixed by the lattice symmetries. A derivation 
of equation (17) can be found e.g in Ref. [52]. We offer 
a condensed review of the derivation in Appendix A. 

In many situations, the factor system Eq. (18) is trivial. 
This is the case, especially, 

(i) for symmorphic lattices because all tj are zero, 

(ii) for momenta k inside the Brillouin zone because 
all gi are zero. 

In these situations, the definition (17) reduces to that 
of an “ordinary” representation. On the other hand, if 


the factor system is non-trivial, Eq. (17) defines a projec¬ 
tive representation. This situation arises for /e-points on 
the surface of the Brillouin zone of non-symmorphic lat¬ 
tices. Interestingly, irreducible projective representations 
of a group can be higher-dimensional than their ordinary 
counterparts. 

This section is structured in the following way. In sub¬ 
section IV B we explain how the symmetries of the space 
group of the symmetric pyrochlore lattice (# 227) pro¬ 
tect a fourfold degeneracy at the X-point of the Brillouin 
zone and why the spectrum disperses linearly around it, 
leading to a Dirac node. In subsection IV C we show how 
breaking inversion symmetry splits the fourfold degener¬ 
acy into two linearly dispersing twofold degeneracies at 
different energies. The lower band of the upper represen¬ 
tation crosses the upper band of the lower representation 
along a 2D sheet in momentum space. In subsection IV D 
we explain why this crossing gaps out everywhere except 
of four points where it is protected by lattice symmetries. 
These four points are the Weyl nodes of the WSM phase 
in our model. 


B. The X point in the symmetric pyrochlore lattice 

1. Summary 

We will now analyze in more detail how the symme¬ 
tries of the non-symmorphic pyrochlore lattice protect 
the Dirac nodes at the X point, where the non-trivial 

application of Eq. (17) arises. Taking the spin-orbit cou- 

_x 

pling into account, the little co-group G contains 32 el¬ 
ements and is isomorphic to the double-valued Ac¬ 
cording to the tables in Ref. [52] , the only irreducible pro¬ 
jective representation compatible with the factor system 
is four-dimensional. This already points to the four-fold 
degeneracy observed in Fig. 2(a). It also implies that at 
a commensurate filling with 2 + 4 n (n E No) electrons 
per unit cell, a band insulator is not possible. This re¬ 
sult complements a similar result found in Ref. [43] for 
the case of a vanishing spin-orbit coupling but with ar¬ 
bitrary electron-electron interaction. 

_x 

By considering the generators of G , it is possible to 
find the factor system and a symmetry-adapted basis that 
spans the 4DIR. This is achieved in Eq. (21). We further 
construct the effective low-energy Hamiltonian, given by 
Eq. (26), which demonstrates that the spectrum indeed 
disperses linearly around the 4DIR. 


2. The factor system 

To be specific, we consider in the following the X point 

_x 

with coordinates X = ^-(1, 0, 0). The little co-group G 
is generated by point operations X, C and C± x where 
+ is an n-fold rotation around axis i in the positive 
(counter-clockwise) direction and the centre of symmetry 


These symmetry operations commute, i.e. 




X 

C 2 + , 

Ctx 



(to) 

(0) 

(to) 

X 

0 9 ) 

-1 

+1 

-1 

C a + 

( 9 ) 

-1 

+1 

-1 


(0) 

+1 

+1 

+ 1 


TABLE I. Factor system of the generators of the little co¬ 
group G . The vector indicated next to each point operation 
in the column is the non-Bravais lattice shift associated with 
it. The vector indicated next to it in the row gives the recip¬ 
rocal lattice shift of the X point under the point operation. 

is the centre of any (but fixed) tetrahedron. Of these 
generators, 

• X and C± x are associated with a non-Bravais lattice 
shift to = —a( 1,1, 1 ), 

• X and C~ 2 Z transform the X point to an equivalent 
point displaced by g = — ^( 1 , 0 , 0 ). 

With this information and knowing that exp (—i g • to) = 
— 1 , it is easy to find the factor system between the group 
generators listed in Table I. 

An important consequence of the non-trivial factor sys¬ 
tem in Table I is that certain commuting point-symmetry 
operations are represented by anXcommuting operators 
and vice versa. We illustrate this fact with the relation 
between space inversion X and the two-fold rotation C^ z . 



FIG. 5. (Color online) The action of {X | to}, where X is the 
inversion around the centre of symmetry indicated by the star, 
and to is a non-Bravais vector shift indicated by the blue ar¬ 
row in the upper left tetrahedron. This symmetry operation 
maps the primed sites to the non-primed ones and is equiv¬ 
alent to the space inversion around site 1. Note that sites 1 
and 2 do not change their x-coordinate under the transfor¬ 
mation while sites 3 and 4 are shifted by —2a. This means 
that {X | to} preserves the amplitude of Bloch wave functions 
with k — X = ^(1,0,0) on sites 1,2, but it changes the 
amplitude on sites 3,4 by a factor exp(i^ • 2a) = — 1. The 
non-uniformity of these factors is a consequence of the non- 
symmorphicity of the symmetry operation. 


loC+=C+ol. (19a) 

However, from Table I it follows that their (projective) 
representations anticommute : 

D X (1)D X (C+) = -D X (C+)D X (1). (19b) 

The origin of Eq. (19b) can be traced back to the non- 
symmorphicity of the pyrochlore lattice. As shown in 
Fig. 5, inversion {X, to} maps three out of four sites to 
neighboring unit cells. As a result, when inversion acts on 
Bloch wave functions at X, non-uniform exponential fac¬ 
tors exp (iX • r ) have to be taken into account. Indeed, 
in the convention of Fig. 5, the amplitudes at sites 1,2 
are unchanged while the amplitudes at sites 3,4 acquire 
a factor 

exp (iX • 2613 ) = —1 = exp (iX • 26 i 4 ). (20) 

The operation C^ z exchanges sites 1,2 with sites 3,4 i.e. 
those that change sign under Dx(X) with those that do 
not. This is just the statement of Eq. (19b). 

3. Construction of a symmetry adapted basis 

We can use our formalism to deduce the existence of 
the fourfold degeneracy at the X point by constructing 
a symmetry-adapted basis. Note that the representation 
D X (C+) = Dx{C^ x ) 2 commutes with the representation 
Dx{X) and that they both commute with XL(X). Hence, a 
basis can be found that diagonalizes these three operators 
simultaneously. Because DxiX) 2 = X^C ^) 2 = — 1 (see 
Table I), both Dx{X) and X^C^) have eigenvalues =bi. 
Let 

IVh) ( 21 a) 

be an eigenvector of 'H(X) with eigenvalue A — ±i under 
Dx{X) and eigenvalue B = d=i under Dx(C 2 },). Then, 
the states 

1 ^ 2 ) =D x {Ci z ) (21b) 

IV> 3 > =D x (Ct x ) |V>i) ( 21 c) 

IV>4> =D X (C+)D X (C+) |^) ( 21 d) 

are eigenvectors at the same energy with distinct eigen¬ 
values under Dx{X) and DxiC^x), as indicated in Ta¬ 
ble II. Hence, the four states (21) are mutually orthogo¬ 
nal and span the (projective) 4DIR at the X point. We 
also see that the 4DIR splits, if any of the symmetry el¬ 
ements {I | to}, {C^ x | to} and { C^ z | 0} is removed from 
the space group. This agrees with findings of Ref. [53] 
that a rotation symmetry is an essential ingredient to ob¬ 
tain a protected Dirac node at a time-reversal invariant 
momentum. Note that the presence of time-reversal T is 
not relevant for the existence of the 4DIR. 
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M> 

M) 

M> 

M) 

Dx(T) 

+A 

-A 

-A 

+A 

Dx(C+) 


-B 

+B 

—B 


TABLE II. Eigenvalues of the states (21) under Dx(X) and 
Dx(C 2 x )• A and B can both assume values ±i. 


4 . Linear dispersion 


We further derive the linear dispersion of the spec¬ 
trum around the X point. A formal group theoretical 
argument, also given in Ref. [42], is expounded in Ap¬ 
pendix B. Here, we instead present a reasoning that ex¬ 
plicitly shows the role of symmetries. 

Insight can be gained by expanding the Hamiltonian 
perturbatively in the momentum p around the X-point. 
The first order term is 


[ftpert..(P)]y = Ml 


dk 


x 


■p) M> 


^ i \[H(X + p)-n{X)]\^ j ). ( 22 ) 


If there is a direction of p for which some of the ma¬ 
trix elements are non-zero, the spectrum might disperse 
linearly in that direction. Otherwise, the spectrum dis¬ 
perses at least quadratically in any direction. 

Symmetries pose constraints on the matrix elements of 
Eq. (22). To find them, we use the transformation laws 
of the Bloch functions (21) deduced from Tables I and II 
together with the fact that the perturbation Hamiltonian 
transforms according to a vector representation, 


R : 


f dRjk) 
\ dk 



i —y 



(23) 


_X 

under a point operation R G G . Equation (23) follows 
from R : H(X + p) •-> H{RX + Rp) = R(X + Rp). 

We start by considering X, which flips the sign of 
all components of the vector p. If the Bloch functions 
{ipi | and | i/jj) have opposite eigenvalues of Tx(I), +A 
and —A, the corresponding matrix element (22) is trans¬ 
formed to minus itself under X and hence must be zero. 
This reasoning forces the matrix elements indicated in 
Table III by the crossed font “p( n to vanish. 

We further consider the operation C^ x which flips the 
sign of p y and p z and preserves the sign of p x . If the 
Bloch functions | and | i/jj) have the same eigenvalue 
±B under C^, they produce a factor (±T>) 2 = — 1 under 
that transformation. This means that the corresponding 
p x matrix element maps to minus itself and must vanish. 
On the other hand, if the two Bloch functions have op¬ 
posite eigenvalues, they produce a factor of +1 and the 
corresponding p y and p z matrix elements are forced to 
be zero. Both of these constraints are indicated by the 
back-crossed terms in Table III. 

The remaining generators of the little co-group do not 
force any of the remaining matrix elements to be zero. 
Since there are p X: p y and p z terms uncrossed for some 



M> 

(+^4, +B) 

M) 

IV’s) 

(—A, +B) 

M) 

(+A,-B) 

Ml (-A-B) 

X,W,X 

>K’Py’Pz 

Px,X’X 

X> XX 

Ml (+A+B) 

>K’Py’Pz 

XX, X 

XX-X 

PxiPy^if^ 

Ml (+A-B) 

Px 1 'Py^ , 

X> XX 

X,X>X 

X’Pv’Pz 

Ml (~A,+B) 

X’XX 


Py 5 Pz 

X’X’X 


TABLE III. The two terms in brackets indicate the eigen¬ 
values of the states under Dx(X) and Dx(C'^ x ), respectively. 
The crossed ” terms [shown in blue and green in the online 
version] indicate matrix elements (22) that vanish due to X 
and the back-crossed terms “pff [shown in blue and green in 
the online version] indicate those that vanish due C} x symme¬ 
try. Since all p x ,p y and p z remain uncrossed for some pair of 
states, the spectrum disperses linearly in all directions around 
the X point. 


pair of wave functions in Table III, the spectrum disperses 
linearly in all directions. 


5. Effective Dirac Hamiltonian 


To deduce the form of the effective Hamiltonian, we 
also analyze how symmetries relate the non-vanishing 
matrix elements in Eq. (22). First, the rotation C^ z leads 
to 

[ftpert iPx )] 13 = ~ [Hpert.(Px )] 24 (24a) 

[ftpert .(Py)\ 12 = + ['Hpert .(P »)] 21 ( 24b ) 

[^pert,(P *)] 12 = - [ftpert.(P *)] 21 ( 24c ) 

and the same relations with (1,2,3,4) «->> (3,4,1,2). Sec¬ 
ond, the remaining little co-group generator C± x , which 


maps (p x ,p y ,p z ) 

1 to (p x , 

Pz , Py) 

, leads 

to 


[^pert. 

. (Py)] 12 

= {-B) 

[Wpert. 

■ (P 2 )] 34 

(25a) 

[Wpert. 

•(P*)]l 2 

= {+B) 

^pert. 

(Py )] 34 

(25b) 

[^pert. 

■(P*)] 13 

= (+B) 

[^pert. 

, (Pa,)] 3 i 

(25c) 

and opposite sign relations with 


M) ** 

(4,3,2,!). 


Combining these results with hermiticity and time- 
reversal leads to an effective Dirac Hamiltonian 


^Dirac(P) = £( X )1 + a Px? 1 + b (p y T 2 + p z V 3 ) (26) 

with a, b G R. The Dirac matrices T^ in our basis are 
given by 


ri = -L (a x + i Ba y ) 0 a z 

(27a) 

r 2 = 1 0 cr x 

(27b) 

T 3 = \B 1 ® (Jy 

(27c) 

r 4 = a z 0 a z 

(27d) 

anticommutation relation 


{ r i,Tj} = 2 Sij. 

(28) 
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The first set of the Pauli matrices in Eq. (27) operates 
on the 2x2 blocks of the matrix Hamiltonian (26), and 
the second set acts within these blocks. Diagonalizing 
the Dirac Hamiltonian (26) leads to a spectrum with two 
doubly degenerate linearly dispersing bands 

e(X + p)= e(X) ± + & 2 (p2 +p2). (29) 


C. The X point in the breathing pyrochlore lattice 

Breaking of the inversion symmetry decreases the al¬ 
lowed degeneracy at the X point from 4 to 2, i.e., the 
4DIR splits into two 2DIRs with the chemical potential 
in-between. The dispersion around each of these nodes is 
linear within the square face of the BZ and quadratic in 
the perpendicular direction. As discussed in Sec. IV D, 
the Weyl nodes appear as symmetry protected crossings 
of the up-dispersing band of the lower 2DIR with the 
down-dispersing band of the upper 2DIR. 

More specifically, the little co-group at the X point 
of the Z-broken pyrochlore lattice is isomorphic to the 
double-valued D^d crystallographic point group. It is 
generated by the rotation C^ z and the improper rotation 

S 4x = I ° c tx '■ ( x > y , z ) ^ (-*> - y ), (30) 

that satisfies (-S^) 2 = C^. Note that the Z-broken py¬ 
rochlore lattice is symmorphic, so the factor system is 
trivial and the degeneracy is given by the “ordinary” IRs. 
According to [52] , all IRs that produce a minus sign un¬ 
der a 27r-rotation are two-dimensional. 

Let us now construct the states that span the 2DIR. 
Note that it is possible to simultaneously diagonal¬ 
ize the Hamiltonian H(X), and operators Dx{S^ x ) and 
D X (C+). Let 

1 4>t) ( 31a ) 

be an eigenstate of H(X) with energy e(X) ± A/2 and an 
eigenstate of DxiC^x) with eigenvalue B = =bi. It follows 
that 

\4)=D x (C+)\^) (31b) 

has the same energy as |0i) but the opposite eigenvalue of 
DxiC^x)- This implies that the two states are orthogonal 
and span the 2DIR at the X point. If an appropriate 
phase of |Vq) in (21a) is adopted, then 



I0i> 

102) 


i+B) 

i-B) 

(0i 1 i-B) 

'ZV’K’K' 

P^,Py,Pz 

(02 1 (+B) 

P^,Py,Pz 



TABLE IV. The ±L> terms indicate the eigenvalue of the state 
under a rotation, which is the square of the eigenvalue un¬ 
der the improper Si, rotation. The blue crossed “jV” terms 
indicate matrix elements (22) that vanish due to C^ x and the 
orange back-crossed terms indicate those that vanish 

due to S± x . The uncrossed black terms are allowed by sym¬ 
metry. The absence of p x implies that the spectrum disperses 
quadratically in this direction. 


where T^ = — ^ [T*, Tj]. The effect of such a term on the 
Dirac Hamiltonian has been discussed in Refs. [20, 36] 
and agrees with our observations, e.g. in Fig. 2(b) and 
Eq. (36) below. 

In analogy with the previous subsection, we can also 
investigate the dispersion of the 2DIRs around the X 
point using the symmetries. They again lead to cer¬ 
tain constraints on the matrix elements between \<pi) and 
| 02 ) : If ((/>i | and | (j>j) have the same eigenvalue un¬ 
der Dx(C f f^), the corresponding p x matrix element must 
vanish. If their eigenvalues are opposite , the p y and p z 
elements are forced to be zero. These matrix elements 
are indicated in Table IV by the crossed “pf” terms. On 
the other hand, the improper rotation sr* flips the sign 
of p x , hence the matrix elements ( 4>i\p x \ fa) are mapped 
to minus themselves under the transformation and have 
to vanish. They are indicated in Table IV by the back- 
crossed terms. 

The remaining matrix elements are non-zero. The im¬ 
proper rotation S± x relates 


^pert .iPy) 
,X,± 


J 12 


= (-B) 


^pert .(.Py) = (+£) Kert .(Pz) 


J 21 


UpeftiPz) 

X,± 


J 12 


J 21 


and the rotation C£ z further leads to 


(34a) 

(34b) 


Kett.iPy) 

- + 

12 

^pei iiPy) 

Kett.iPz) 

12 

Kett.iPz) 


(35a) 

(35b) 


kb2> = ^(lV’i,2>±i’/B|V’3,4>) (32) 

where +(—) sign refers to the upper (lower) 2DIR, as¬ 
suming that A > 0, and the branch cut of the square 
root is along the negative real axis. As a consequence, in 
linear order in A and p, the Dirac Hamiltonian (26) ac¬ 
quires an additional term due to the inversion-symmetry 
breaking, 

^Dirac (p) ^ ^Dirac (p) + ^ F 14, (33) 


Time reversal does not lead to further constraints. Af¬ 
ter considering the hermiticity, the perturbation Hamil¬ 
tonian can be written as 

iP) = ( £ ( x ) ± f) 1 + b iPyVx + iBp z a y ). (36) 

with c G R. The two bands disperse linearly along the 
p y and p z directions and quadratically along the p x di¬ 
rection. This is in agreement with analysis of Eq. (33) 
in Ref. [36]. 
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The breaking of the inversion symmetry splits the lin¬ 
early dispersing 4DIR discussed in subsection IV B into 
the two 2DIRs discussed in IV C that disperse linearly 
only within the square face of the Brillouin zone and 
quadratically along the TX direction. The up-dispersing 
band of the lower 2DIR crosses the down-dispersing band 
of the upper 2DIR on a 2D sheet that locally looks like a 
one-sheet hyperboloid with the X point in the centre. At 
a general k- point, this crossing gaps out because there 
is no symmetry to protect it. Potential exceptions are 
k -points with higher symmetry and we therefore discuss 
the XU and the XW lines in more details. 


1. The XU-line 

Upon breaking the inversion symmetry, the four bands 
originating from the 4DIR at the X point disperse as 
indicated in Fig. 6(a). The question is whether the band 

crossing is protected or gapped. 

_xu 

The little co-group G is cyclic, the only generator 
is a mirror symmetry The possible eigenvalues of 

DxuC^di) are ±i and correspond to two different IDIRs. 
Two bands with the same eigenvalue can hybridize and 
their crossings generally gap out. On the other hand, 
two bands with opposite eigenvalues cannot hybridize and 
hence their crossings are protected. To ascribe the ap¬ 
propriate 1DIR to each of the four bands illustrated in 
Fig. 6(a), we devise the following two arguments. 

First, if we restore the inversion symmetry, then T o 
X enters the little co-group, and the bands 1 and 2 in 
Fig. 6(a) become degenerate. This element commutes 
with (Jdi- Furthermore, both Tol and cr^i are associated 
with zero shift g = 0 in momentum space, leading to a 
trivial factor system, therefore 

[E>xu{T oX),D X u(<Xn)] = 0. (37) 

Note that Dx\j(T ol) is antiunitary because of the time- 
reversal. Now, if state |?/q) from band 1 is eigenvector 
of Dyjj(<Jdi) with eigenvalue B = ±i, then the state 
|^ 2 ) = E>xu{T o T) |^i) from band 2 has eigenvalue 
B* = —B under Dxu( cr di) and hence belongs to the 
other representation. Similar conclusion can be found 
for bands 3 and 4. 

Second, if we keep inversion symmetry broken and 
move along the XU line to the X point, a perpendicu¬ 
lar mirror symmetry 0^2 enters the little co-group that 
fulfils Gd\ 0 <?d 2 = E o (j d 2 0 &di ? where E is a 27r-rotation. 
Therefore 

{Dx{vdi),Dx{cFd2)} = 0. (38) 

The state 1 ^ 3 ) = Dx{cF d 2 ) |Vh) from band 3 has the same 
energy as |^i) but has eigenvalue —B under Dx(<Jdi)- 
This means that bands 1 and 3 belong to different IDIRs. 
The same conclusion can be found for bands 2 and 4. 




(a) (b) 

FIG. 6. (Color online) (a) A schematic sketch of the four 
bands along the XU lines, originating from the 4DIR at the 
X points. There are only two IDIRs along this line. Bands 
belonging to one of them are indicated by dashed red lines, 
those belonging to the other representation are shown in solid 
blue. The bands that cross belong to the same representation, 
hence the crossing gaps out. (b) Analogous analysis along the 
XW line reveals that the crossing bands belong to different 
representations. Such a crossing is protected by symmetry 
and yields a Weyl node that can be gapped out only by an¬ 
nihilation with another Weyl node. 

Our conclusions are summarized in Fig. 6(a). Bands 
with eigenvalue +B are indicated by dashed red lines, 
while those with eigenvalue —B by solid blue lines. We 
see that the two crossing bands belong to the same rep¬ 
resentation, so the crossing will in general be gapped out 
by hybridization. 


2. The XW-line 

We finally investigate the crossing along the XW line. 

For concreteness, we consider the line parallel with the 

—xw 

z-axis. The little co-group G is generated by a two¬ 
fold rotation C^ z - Possible eigenvalues of DxwiC^) are 
±i and correspond to two different IDIRs. We want to 
assign a representation to each of the bands sketched 
in Fig. 6(b). The argumentation proceeds again in two 
steps. 

First, if we restore the inversion symmetry X, then ToX 
and mirror symmetries a x and a y appear in the little co¬ 
group, and bands 1 and 2 become degenerate. Further¬ 
more, the system is again non-symmorphic and the factor 
system is non-trivial. A careful calculation reveals that 
representations Dxw(C^), ftw(^) and D X wK) mu¬ 
tually commute, so a basis can be found that diagonalizes 
all of them simultaneously. 

Let l^i) be an element of such a basis with eigenvectors 
B = =bi, A x = ±1, A y = =bi, respectively. Since 

^xw(^)^xwK) -Pxw (C^) (39) 

the eigenvalues are constrained by A x A y = — B. It turns 
out that the state ^ 2 ) = Dxw{T ol) |?/q) from band 
2 has the same energy and eigenvalues B,—A x ,—A y , 
respectively. The unchanged sign of the eigenvalue of 
Exw{G 2 Z ) means that bands 1 and 2 belong to the same 
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representation. Analogous statement can be made about 
bands 3 and 4. 

Second, if we move towards the X point, then 0% en¬ 
ters the little co-group and bands 1 and 3 become degen¬ 
erate. Since C^. o C^y = E o C^y ° C^, we find that 

{D X (C+),D X (C+)} = 0. (40) 

This means that the state 1^3) = Dxi/C^y) from band 3 is 
an eigenstate at X with the same energy and eigenvalue 
—B under Dx{C 2 Z ) This implies that bands 1 and 3 in 
Fig. 6(b) belong to different representations. Analogous 
conclusion can be derived for bands 2 and 4. 

Note that the bands crossing along the XW line belong 
to different representations. This means that the crossing 
is protected. It corresponds to one of the Weyl nodes of 
the Weyl semimetal phase in our model. The other Weyl 
points are related by crystal symmetries. 

V. {111} AND {111} SURFACE STATES 
A. Terminations of a (111) slab 

The topological nature of the WSM phase is mani¬ 
fest in the exotic surface states that have the form of 
Fermi arcs and that are robust against all local pertur¬ 
bations [15, 17]. The endpoints of the Fermi arcs are 
given by projections of the bulk Weyl nodes onto the 
surface Brillouin zone (SBZ). In the following, we focus 
on the {111} and {111} surfaces as these are the natu¬ 
ral cleavage planes of the pyrochlore oxides. Besides the 
topological surface states exhibited by the WSM of model 
(12), we also identify non-topological surface states in the 
DSM and the INS phase for certain terminations of the 
pyrochlore lattice. 

Along the {111} directions, the pyrochlore lattice can 
be viewed as a stack of alternating layers of triangu¬ 
lar and kagome lattices. Hence, the two simplest open 
boundary conditions correspond to terminating the crys¬ 
tal at sites making up either the kagome lattice (K) or the 
triangular lattice (T). Note that the kagome termination 
cuts the outermost tetrahedra. Because of the inequiva¬ 
lent bonds in the Z-broken state, we further have to spec¬ 
ify whether the bonds connecting the outermost layer of 
sites to the next layer are strong (S) or weak (W). This 
gives four possible terminations illustrated in Fig. 7. For 
the DSM phase the labels S and W are redundant and 
can be dropped. A suitable geometry for experimental 
studies is a thin slab with large lateral dimensions [23]. 
If we set the (111) surface on the top and we limit our 
attention to (5 > 0, then only KS and TW can be realized 
on the top, and only KW and TS on the bottom surface, 
as shown in Fig. 7. 

The SBZ has the shape of a regular hexagon and its 
construction is indicated in Fig. 8(a). Since the twelve 
Weyl nodes project onto separate points, six distinct 
Fermi arcs are expected. Upon increasing the staggered 



FIG. 7. (Color online) Illustration of the four crystal termi¬ 
nation in the {111} and {111} directions considered in Sec. V. 
The crystal can terminate either at the triangular (T) or the 
kagome (K) lattice, and for S ^ 0 we further have to specify 
whether the bonds connecting the outermost layer of sites to 
the next one are strong (S) or weak (W). Only TW and KS 
terminations are possible on the top, and only TS and KW 
are possible on the bottom of a sample in the slab geometry. 


strain S, the projections of the Weyl nodes move to the in¬ 
side of the SBZ along the trajectories shown in Fig. 8(b). 
A summary of the surface states for the four different 
terminations shown in Fig. 7 and varying values of S [the 
remaining parameters are fixed as in Eq. (13)] is provided 
in Fig. 9, where we show density plots of the the surface 
spectral function. 


B. Strain dependent surface states 

1. Surface states of the Dirac semimetal 

The column 5 = 0 in Fig. 9 corresponds to the DSM 
phase. We find no surface states for the triangular ter¬ 
mination. For the kagome termination, we find a single 
non-degenerate Fermi line connecting the M points of 
the hexagonal SBZ, which are the projections of the bulk 
Dirac nodes. 

Since the surface breaks symmetries that protect the 
bulk Dirac nodes, a small gap, which shrinks with in¬ 
creasing system width, opens in the (111) film. As a 
consequence, the surface states for the kagome termina¬ 
tion of the DSM phase do not exactly pass through the 
M points but slightly avoid them. In general, for band 
insulators, the Fermi line has to encircle an even or odd 
number of time-reversal invariant momenta (TRIMs) of 
the SBZ, depending on the strong Z 2 invariant of the 
band structure [54, 55]. We therefore expect that either 
the situation shown in Fig. 8(c) (normal insulator) or in 
Fig. 8(d) (topological insulator) is realized. We numer¬ 
ically find that the first possibility (normal insulator) is 
realized, which is consistent with the fact that there are 
no surface states for the triangular termination. 
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FIG. 8. (Color online) (a) Construction of the surface Bril- 
louin zone (SBZ) along the {111} directions. The time- 
reversal invariant momenta (TRIMs) of the SBZ are given 
by projections of the bulk TRIMs. The projections of the 
Weyl nodes of the two possible chiralites are also indicated, 
(b) Definition of the high-symmetry points T, M and K in 
the SBZ. The three inequivalent M points and the T point 
are TRIMs. The solid red and dashed green lines indicate the 
trajectory of the projections of the Weyl nodes of opposite 
chiralities upon increasing the value of 5. The crosses indi¬ 
cate where the Weyl points annihilate for the critical 6 C . (c) 
and (d) The DSM phase gets gapped for finite thickness of 
the system. The Z 2 invariant of the insulating phase depends 
on the number of TRIMs encircled by the surface Fermi lines. 
Situation (c) corresponds to a normal and (d) to a topological 
insulator. 


2. Surface states of the Weyl semimetal 

In the WSM phase, the surface states form Fermi arcs 
connecting the projections of the Weyl nodes. Interest¬ 
ingly, the connectivity of the Fermi arcs, i.e. the way the 
Weyl nodes of opposite chirality are paired into the arcs, 
depends on the termination of the sample. In Ref. [56], a 
similar dependence appears naturally in a toy model con¬ 
sisting of a stack of alternating electron and hole Fermi 
surfaces. Here, we directly observe this phenomena in a 
microscopic model. 

For the TW termination, we observe six Fermi arcs 
developing between the projections of the Weyl points in 
neighboring SBZs as 5 increases from zero. At the critical 
6 C , three Fermi arcs at the time form a closed Fermi line 
enclosing K and K’, respectively. For even larger values 
of 6, all surface bands are shifted away from the bulk 
chemical potential. 


For the KS termination, the very opposite happens. 
The closed Fermi lines present in the DSM phase splits 
into six distinct Fermi arcs upon breaking the inversion 
symmetry. Increasing the staggered strain S leads to a 
shrinking of the Fermi arcs and to a complete disappear¬ 
ance of the surface states at S c . 

For the TS termination and S > 0.04, the same pairs 
Weyl nodes as for the KS termination are connected but 
with opposite curvature. They also disappear at S c . In 
Subsec. V C, we separately discuss the interesting param¬ 
eter range 0 < 5 < 0.04. 

Finally, for the KW termination, additional Fermi arcs 
appear very close to the Fermi lines present in the DSM 
phase. This fact makes their observation obscured in 
Fig. 9. For the INS phase we find a pair of closed Fermi 
lines encircling the T point. 

3. Surface states of the insulator 

In the INS phase, we observe no surface states for the 
TS, TW and KS terminations. However, for the KW ter¬ 
mination, there are two surface bands crossing the bulk 
chemical potential that encircle the T point. These ob¬ 
servations are in accordance with the trivial Z 2 invariant 
of the INS phase. 

One possibility to understand the “non-topological” 
surface states of the KW termination is to consider the 
(unphysical) limit 5 1 and p 0 which corresponds 

to a lattice of isolated tetrahedra. A lattice of isolated 
tetrahedra exhibits only flat bands. The boundary con¬ 
sists of 

(i) isolated tetrahedra for the TS and KS terminations, 

(ii) isolated triangles for the KW termination, 

(iii) isolated points for the TW termination. 

In the situation (i), the bulk and the surface have iden¬ 
tical spectra so the chemical potential lies in the gap of 
both. For sufficiently weak coupling between the tetrahe¬ 
dra, this observation has to remain valid and we do not 
expect surface states for the TS and KS terminations. 

For the other two situations, we need to know the 
eigenenergies of the corresponding boundary objects. Fo¬ 
cusing on s = +1, the states of an isolated tetrahedron 
with spin-orbit coupling parametrized by R lie at energies 

4 4) «6, 4 4) = -2(l —4 R), e$ = -2(l + 2i?). (41) 

Due to time-reversal symmetry, all levels are doubly de¬ 
generate. It follows that for R = —0.4 (as in Fig. 9), the 
chemical potential satisfies e^ 4 ] < fi < e^. On the other 
hand, the energy levels of an isolated triangle appearing 
at the KW termination are 

4 3) = 1 + 2 R± x/3 [3 + 4(i? — l)i?], (42) 

4 3) = -2(l + 2 R). (43) 
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FIG. 9. (Color online) Surface spectral function within the SBZ defined in Fig. 8(a), plotted for varying value of 5 and varying 
surface terminations illustrated in Fig. 7. The blue circles in the DSM column represent projections of the bulk Dirac nodes 
into the SBZ. In the same way, the red square and the green triangle markers in the WSM columns represent the Weyl nodes of 
the two opposite chiralities, and the pale crosses in the WSM-to-INS column represent the fc-points where pairs of Weyl nodes 
annihilate. In the INS phase, the chemical potential is set to the middle of the bulk gap. 


Hence, the highest energy level also satisfies < 

which pins fi at for T = 0 in the charge 
neutral system. The eigenenergy of an isolated point 
appearing at the TW termination is simply 

4 1} = 0, (44) 

which too lies between e^ 4 ] and e 1 , 4 " 1 for R = —0.4. We 
therefore expect that there are surface states in the bulk 
gap both for the KW and the TW termination, even if 
we reintroduce the coupling between the tetrahedra. We 
indeed do observe this. In Fig. 9, we fixed the chemical 
potential in the INS phase at the bulk value. For the KW 
termination, the surface states cross the bulk chemical 
potential, for the TW termination they lie below it. 


C. Weyl-Lifshitz transitions 

As illustrated in Fig. 9, the connectivity and shape 
of the Fermi arcs in the WSM phase depend on the 
termination of the sample. For the studied surfaces of 
our model, there are three possible connectivities that 
respect the symmetries of the system and that do not 


contain Fermi arcs crossings (which are expected to be 
gapped out) [Fig. 10(a-c)]. Interestingly, it is possible 
to continuously move from one connectivity to another 
one, giving rise to a Lifshitz transition of the Fermi arcs 
that we dub Weyl-Lifshitz transition. More precisely, the 
Weyl-Lifishitz transition is characterized by a topological 
change of the closed Fermi lines formed by the Fermi arcs 
of the top and bottom surfaces. 

We have observed such a transition for the TS ter¬ 
mination for varying staggered strain: For small S, the 
connectivity of the Fermi arcs is identical to that of the 
TW termination, see Fig. 9. At S « 0.015, a part of 
the topological surface band crosses the chemical poten¬ 
tial from below at the two K points. The corresponding 
Fermi lines grow and at S « 0.04 they touch the origi¬ 
nal Fermi arcs. A reconnection of the Fermi arcs from 
situation in Fig. 10(b) to that of Fig. 10(a) occurs. 

Another possibility is to tune the shape and connectiv¬ 
ity of the Fermi arcs by applying a surface gate potential. 
We model such an experiment in a simplified manner 
by adding an on-site potential V to the outermost layer 
of atoms. This effectively shifts the energy of the sur¬ 
face states and those that were originally away from the 
Fermi level can be tuned to cross it. The surface states 
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are then allowed to hybridize with the Fermi arcs, lead¬ 
ing to a reconnection of the Weyl nodes. Figure 11(a) 
shows the transition between the connectivity shown in 
Fig. 10(a) and the one shown in Fig. 10(b) by applying 
a positive surface potential to the KS termination. Sim¬ 
ilarly, we observe a transition from the situation shown 
in Fig. 10(a) to the one shown in Fig. 10(c) if we apply 
a negative surface potential. 

A fascinating aspect of the Weyl-Lifshitz transition is 
that it changes the number of times an electron crosses 
the bulk in order to complete its semiclassical orbit in an 
external magnetic field. For example, if we assume the 
TS termination on the bottom and the KS termination on 
the top surfaces, the transitions shown in Fig. 11 changes 
the number of bulk crossings from 2 to 6 in (a) and from 
2 to 12 in (b). 

Such a dramatic change might be observable in quan¬ 
tum oscillation experiments in very clean samples for 
which the mean-free path l exceeds the system width 
L. In fact, according to Ref. [23], the surface-state re¬ 
sponse shows periodic-in- 1/B oscillations for fixed /i in 
the resistivity (Schubnikov-de Hass effect) or magneti¬ 
zation (de Haas-van Alphen effect), which results from 
energy levels periodically crossing the chemical potential 
[i. In the semiclassical approximation, this happens for 




FIG. 11. (Color online) Weyl-Lifshitz transitions at 6 = 0.08 
realized by tuning the on-site potential on the outermost layer 
of atoms for the KS termination. Assuming the TS termina¬ 
tion on the opposite surface according to Fig. 10(d), we find 
that (a) the transition at positive V = 0.93 changes the num¬ 
ber of times the electron orbits cross the bulk from 2 to 6, and 
(b) the transition at negative V = —0.82 changes the same 
number from 2 to 12. 


the nth energy level if 


1 


e7rv(n + 7 ) 
ko/i 


I COS (j) | 


ecL 
hko ' 


(45) 


Here, fco denotes the Fermi arc length, v the characteris¬ 
tic Fermi velocity, the angle of the magnetic field with 
respect to the surface normal and L the thickness of the 
sample. 7 is of order unity and incorporates small n 
quantum effects and 2 c denotes the number of bulk cross¬ 
ings (in Ref. [23], only c— 1 has been considered). Hence, 
by measuring the dependence on the field-direction, one 
can extract the second term of Eq. (45). In an ex¬ 
periment with a surface gate potential similar to Fig. 11, 
the change in originates from a change of the length 
of the Fermi arcs Ako but also from the change in the 
number of bulk crossings Ac. 


FIG. 10. (Color online) (a,b,c) Three connectivities of the 
Fermi arcs that preserve the symmetry of the system. The 
dark red and pale green points correspond to the projections 
of the bulk Weyl nodes of opposite chiralities. The arrows 
indicate direction of motion of the electrons in the quantum 
oscillations experiment proposed in [23] when a magnetic field 
is applied perpendicular to the surface, (d) The oriented ver¬ 
tical lines correspond to the chiral Landau levels emanating 
from the Weyl nodes, and act as “conveyor belts” transport¬ 
ing the electrons between the top and bottom surfaces of the 
slab. If the indicated connectivity of the Fermi arcs is real¬ 
ized on the bottom surface, the electron orbits traverse the 
bulk 2, 6 and 12 times, respectively, if connectivities (a),(b) 
and (c) are realized on the top surface. Realizing situations 
(b) on the top and (c) on the bottom leads to electron orbits 
traversing the bulk four times. 


VI. CONCLUSION 

In summary, we studied an effective lattice model of 
spin-orbit coupled electrons on the pyrochlore lattice at 
a commensurate filling that realizes a Dirac semimetal 
phase protected by the non-symmorphic space group. 
Upon coupling the electronic degrees of freedom to the 
lattice, a staggered strain, which breaks the inversion 
symmetry, can spontaneously develop. For increasing 
strain, each Dirac node splits into four Weyl nodes and 
a topological semimetal phase is realized. At a critical 
strain, Weyl nodes of opposite chirality annihilate and 
the system enters a trivial insulating phase. We iden¬ 
tified several interesting features of our model such as 
a reentrant behavior of the Weyl semimetal phase in the 
elasticity vs. temperature phase diagram and the appear- 
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ance of a Griffiths wing in the presence of a staggered 
stress. 

We furthermore presented a detailed group-theoretical 
analysis of the electronic spectrum, which is independent 
of the details of the considered lattice model and only 
relies on symmetries. We highlighted the importance 
of non-symmorphicity for realizing a Dirac semimetal 
and displayed symmetry-based arguments for the appear¬ 
ance of the Weyl nodes in the inversion-symmetry broken 
phase. 

Eventually, we studied the surface states of our model. 
Most notably, we found a dependence of the connectivity 
and the shape of the Fermi arcs on the termination of 
the sample. We also demonstrated that it is possible to 
continuously change the connectivity of the Fermi arcs 
through a Weyl-Lifshitz transition by applying a surface 
potential. A change in the connectivity may also change 
the number of times an electron crosses the bulk in order 
to complete its semiclassical orbit, which leads to a clear 
signal in a quantum oscillation experiment with a varying 
field direction. 

The model considered in this paper is relevant for py- 
rochlore oxides with heavy transition-metal ions such as 
iridates. For these systems, one possibility to tune the 
chemical potential to the Dirac nodes is to consider al¬ 
loys of the type A.2-x^>x^20j. In transition-metal oxides, 
electronic correlations are often important but in 5d sys¬ 
tems, they are in general less pronounced. Moreover, the 
fundamental Berry curvature structure around the Weyl 
nodes is perturbatively stable against interactions [57]. 
In our study, we have therefore neglected these effects and 
instead studied a non-interacting system. However, the 
interplay between electron-electron and electron-lattice 
interactions is a challenging but interesting research di¬ 
rection which we leave for further studies. 

In conclusion, we presented a route in which a non¬ 
magnetic Weyl semimetal is realized as a thermodynamic 
phase with spontaneous inversion symmetry breaking 
upon coupling the electronic degrees of freedom to the 
lattice. Many of the observed phenomena of our model, 
such as the reentrant behavior of the Weyl semimetal, 
the dependence of the surface Fermi arcs on the termi¬ 
nation and the Weyl-Lifshitz transitions are expected to 
be independent of the details of the system and should 
be applicable to a large class of materials, in particular 
including other systems belonging to the #227 (Fd3m) 
space group such as diamond, /3-cristobalite and spinel 
oxides structures. 
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Appendix A: Projective representations of space 
groups 

Here, we provide a derivation of equation (17). A more 
complete discussion can be found in Ref. [52] . For brevity, 
we write IR for irreducible representation and nDIR for 
n-dimensional IR throughout the appendix. 

Every element of a space group G can be expressed 
as a point operation R followed by a uniform shift by a 
vector t which we write compactly as {R | £}, 

{R 1 1} : v Rv + t (Ala) 

Action of the element on functions in real space, including 
Bloch wave functions, is 

{R\t} ■ f(v) h> /[R-'iv-t)} (Alb) 

and the composition rule is 

{i ? 2 | V 2 } O {#! | m = {i ?2 ° -Ri | R2V1 + V 2 }. (Ale) 

The identity element is {E | 0}. 

The set T of all pure translations by Bravais lattice 
vectors is a subset of G, and can be used to write all 
elements of G compactly as 

G = T o {ith | ti} + ... + T o {R n | t n } (A2) 

where the point operations R \,..., R n act at the same 
point in real space and are all different. If the vectors 
ti can all be made zero by a proper choice of the point 
of symmetry, the lattice is called symmorphic. If this 
cannot be done, the lattice is dubbed non-symmorphic. 
Note also that the set 

F = {Ri,..., Rn} (A3) 

of the point operations is always a group, while the coset 
representatives {Ri \ tf} form a group if and only if the 
space group is symmorphic. 

If we take spin-orbit coupling into account, a 2tt- 
rotation of a wave function in real space is accompanied 
by a 27r-rotation of the electron spin which results in 
a sign change of the wave function. This operation is 
not equivalent to identity and we denote it as E. A dvr- 
rotation is equivalent to identity E. 

We adopt the periodic boundary conditions. Then the 
group T is Abelian and as such it has only IDIRs p k la¬ 
belled by momenta k. The corresponding representation 
space is spanned by a Bloch wave function at k 

#a,fe(^) = exp (ifc • r)u a , k (r) (A4) 

where u a , k (r) is the cell-periodic part. The representa¬ 
tions of a pure translation by vector v is 

p k {{E | v}) = exp (-i k • v). (A5) 

The IRs of the space group G can be more than one¬ 
dimensional, but they reduce to the IDIRs (A5) on the 
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subgroup T. Let us consider a representation p of G 
that contains pk in its decomposition on subgroup T, 
i.e. it contains fj a ,k{ r ) as one of the basis vectors in its 
representation space. Then element {R \ t} transforms a 
Bloch function at k into 

{R 1 1 } : exp (ife • r)u a: k{r) 

i-^exp [i k • (R _1 r) — i k • — t) (A6) 

which can be easily recognized as a Bloch function at 
Rk. This means that p necessarily also contains pRk in 
its decomposition on T. 

To find the allowed spectrum degeneracies at k we have 
to consider only those symmetry operations that leave 
the momentum of a Bloch function invariant (modulo 
reciprocal lattice vectors). We construct it as follows. 
Let us denote the subgroup of point operations F that 

_ k 

leave k invariant (called the little co-group of fe) as G . 
Then the group we are looking for is 

G k = {jTo{R i \t l }, Ri£G k . (A7) 

i 

It is a subgroup of G called the little group of k. 

The IRs pk of G k reduce on the subgroup T to IDIRs 
labelled by the same momentum k , so according to equa¬ 
tion (A5) for a Bravais vector v 

pk {{E | v}) = exp (-i k • v)t (A8) 

where 1 is the unit matrix. This means that represen¬ 
tations of Bravais translations commute with representa¬ 
tions of all other elements of G k . 

It is useful to perform a substitution 

p k ({R 1 1}) = exp (-i k • t)D k {{R | £}). (A9) 

The composition rule (Ale) and the representation of 
Bravais translations (A8) imply that 

Dk{{Ri | ti})Dk{{Rj I tj}) 

= exp [—ig { -tj)D k ({R k \ t k }) (A10) 

where g i = ( R~ 1 k) — k is a reciprocal lattice vector, 

_ k 

Rk = Ri o Rj is a point group operation from G , and 
tk is a vector appearing together with Rk in expansion 
(A2). Note that instead of considering the function Dk 
on elements of G fe , we might restrict our attention to its 

_ k 

values on elements of G by defining 

Dk(Ri) •= Dk{{Ri | U}), (Alla) 

Dk(Ri)D k (Rj ) = exp (-igi • tj)D k (Ri o Rj). (Allb) 

This completes the derivation of equation (17). 

The factors exp {—igi • tj ) are completely fixed by the 
crystal symmetry and are referred to as the factor system 
of the projective representation. The projective represen- 

_ k 

tations of group G can be found as ordinary represen¬ 
tation of some larger group that we will refer to as the 


extension group of G . One only has to pick up those 
representations of the extension group that are compat¬ 
ible with the factor system. Reference [52] goes through 
all high symmetry point of all space groups, gives the 
appropriate extension of every little co-group, provides a 
complete list of their representations, and picks up those 
that are compatible with the factor system. 

Appendix B: Linear dispersion around the Dirac 
node 

In subsection IV B we analysed the dispersion around 
the four-fold degeneracy at the X point by studying con¬ 
straints posed on the matrix elements (22) by the gen- 
_ x 

erators of G . In this appendix, we briefly introduce 
an alternative procedure mentioned in Ref. [42]. This 
method is computationally very efficient, but it conceals 
the role of individual symmetries as well as the specific 
form of the k • p expansion (26). 

According to the Clebsh-Gordan decomposition, the 
number of times, a CT , that the complex conjugate repre¬ 
sentation is contained in the product representation 
p M x p u , is given by 

a a = (Bl) 

where C k are the classes of the group G (in our case, 
_ x 

the extension group of G ), \G\ is the number of ele¬ 
ments of G, Nk is the number of elements in C \ c , and 
X stands for the characters of the representations. The 
right-hand side of Eq. (Bl) is symmetric under permuta¬ 
tions of (/i,z/, a), therefore [58] a a is also the number of 
times appears in p v x and the number of times 
p u * appears in p M x p a . 

We further consider the selection rules. Let |^) be dt 1 
states that belong to IR and be d v states that 
belong to IR p v . The scalar products \<l>j) are invari¬ 
ant under G. Using the Schur orthogonality relation, we 
find 

< W> = A X X p*M<r je (g) mrr 

1 1 ki geG 

= j(B2) 

k 

i.e. the products are non-zero iff the functions |^f) and 
1} belong to the same representation. By the symme¬ 
try of (Bl), this is equivalent to stating that the trivial 
representation p 1 is contained in the product representa¬ 
tion p^ x p v . 

We further consider matrix elements 

mo\rp 


(B3) 
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of an operator O that transforms according to represen¬ 
tation p a . It is useful to define states \^ X ^) = Okj \(/>j) 
that transform according to p a x p u . The findings of the 
previous paragraph applied to states |^f) and \^ XV ^) 
imply that matrix elements (B3) are non-zero iff p a is 
contained in the product p M x p v [58]. 

We apply these findings to matrix elements (22). As 
discussed therein, the perturbation Hamiltonian trans¬ 
forms according to the vector representation p vec , and the 
states transform according to representation p 4D . Note 
that we can decompose the product p 4D x p 4D as a sum of 
the symmetric part [p 4D x p 4D ] acting in the hermitian 


sector of the matrix elements, 

Sij (p) = ['Hpert. (p)] ij + [ftpert. (p)] > ( B4 ) 

and the antisymmetric part {p 4D x p 4D } acting in the 
antihermitian sector of the matrix elements, 

An (P) = [H* r (p)].. - (p)] *.. (B5) 

Since, in our case, the two sets of wave functions are 
identical , the antisymmetric sector completely vanishes. 
Hence, the linear correction to the spectrum is non-zero 
iff p vec is contained in the symmetrized Kronecker product 

[p 4D xp 4D ]. 
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